Discovery of novel Mamastroviruses in Bactrian camels and dromedaries reveals complex recombination history.

Virus emergence may occur through interspecies transmission and recombination of viruses coinfecting a host, with potential to pair novel and adaptive gene combinations. Camels are known to harbor diverse ribonucleic acid viruses with zoonotic and epizootic potential. Among them, astroviruses are of particular interest due to their cross-species transmission potential and endemicity in diverse host species, including humans. We conducted a molecular epidemiological survey of astroviruses in dromedaries from Saudi Arabia and Bactrian camels from Inner Mongolia, China. Herein, we deployed a hybrid sequencing approach coupling deep sequencing with rapid amplification of complementary deoxyribonucleic acid ends to characterize two novel Bactrian and eight dromedary camel astroviruses, including both partial and complete genomes. Our reported sequences expand the known diversity of dromedary camel astroviruses, highlighting potential recombination events among the astroviruses of camelids and other host species. In Bactrian camels, we detected partially conserved gene regions bearing resemblance to human astrovirus types 1, 4, and 8 although we were unable to recover complete reading frames from these samples. Continued surveillance of astroviruses in camelids, particularly Bactrian species and associated livestock, is highly recommended to identify patterns of cross-species transmission and to determine any epizootic threats and zoonotic risks posed to humans. Phylogenomic approaches are needed to investigate complex patterns of recombination among the astroviruses and to infer their evolutionary history across diverse host species.

Quasispecies Fitness Partition to Characterize the Molecular Status of a Viral Population. Negative Effect of Early Ribavirin Discontinuation in a Chronically Infected HEV Patient.

The changes occurring in viral quasispecies populations during infection have been monitored using diversity indices, nucleotide diversity, and several other indices to summarize the quasispecies structure in a single value. In this study, we present a method to partition quasispecies haplotypes into four fractions according to their fitness: the master haplotype, rare haplotypes at two levels (those present at <0.1%, and those at 0.1−1%), and a fourth fraction that we term emerging haplotypes, present at frequencies >1%, but less than that of the master haplotype. We propose that by determining the changes occurring in the volume of the four quasispecies fitness fractions together with those of the Hill number profile we will be able to visualize and analyze the molecular changes in the composition of a quasispecies with time. To develop this concept, we used three data sets: a technical clone of the complete SARS-CoV-2 spike gene, a subset of data previously used in a study of rare haplotypes, and data from a clinical follow-up study of a patient chronically infected with HEV and treated with ribavirin. The viral response to ribavirin mutagenic treatment was selection of a rich set of synonymous haplotypes. The mutation spectrum was very complex at the nucleotide level, but at the protein (phenotypic/functional) level the pattern differed, showing a highly prevalent master phenotype. We discuss the putative implications of this observation in relation to mutagenic antiviral treatment.

Deep mutational learning predicts ACE2 binding and antibody escape to combinatorial mutations in the SARS-CoV-2 receptor-binding domain.

The continual evolution of SARS-CoV-2 and the emergence of variants that show resistance to vaccines and neutralizing antibodies threaten to prolong the COVID-19 pandemic. Selection and emergence of SARS-CoV-2 variants are driven in part by mutations within the viral spike protein and in particular the ACE2 receptor-binding domain (RBD), a primary target site for neutralizing antibodies. Here, we develop deep mutational learning (DML), a machine-learning-guided protein engineering technology, which is used to investigate a massive sequence space of combinatorial mutations, representing billions of RBD variants, by accurately predicting their impact on ACE2 binding and antibody escape. A highly diverse landscape of possible SARS-CoV-2 variants is identified that could emerge from a multitude of evolutionary trajectories. DML may be used for predictive profiling on current and prospective variants, including highly mutated variants such as Omicron, thus guiding the development of therapeutic antibody treatments and vaccines for COVID-19.

A key F27I substitution within HCDR1 facilitates the rapid maturation of P2C-1F11-like neutralizing antibodies in a SARS-CoV-2-infected donor.

Although thousands of anti-SARS-CoV-2 monoclonal neutralizing antibodies (nAbs) have been identified and well characterized, some crucial events in the development of these nAbs during viral infection remain unclear. Using deep sequencing, we explore the dynamics of antibody repertoire in a SARS-CoV-2-infected donor, from whom the potent and broad nAb P2C-1F11 (the parent version of Brii-196) was previously isolated. Further analysis shows a rapid clonal expansion of some SARS-CoV-2-specific antibodies in early infection. Longitudinal tracing of P2C-1F11 lineage antibodies reveals that these elite nAbs were rare. Using sequence alignment, structure modeling, and bioactivity analysis based on site-mutated assay, we demonstrate that a key substitution F27I in heavy chain contributes significantly to the maturation of P2C-1F11-like antibodies. Overall, our findings elucidate the developmental process and maturation pathway of P2C-1F11, providing some important information for the design of novel immunogens to elicit more potent nAbs against SARS-CoV-2 infection.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations.

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Comparative analysis of SARS-CoV-2 quasispecies in the upper and lower respiratory tract shows an ongoing evolution in the spike cleavage site.

Studies are needed to better understand the genomic evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study aimed to describe viral quasispecies population of upper and lower respiratory tract by next-generation sequencing in patients admitted to intensive care unit. A deep sequencing of the S gene of SARS-CoV-2 from 109 clinical specimens, sampled from the upper respiratory tract (URT) and lower respiratory tract (LRT) of 77 patients was performed. A higher incidence of non-synonymous mutations and indels was observed in the LRT among minority variants. This might be explained by the ability of the virus to invade cells without interacting with ACE2 (e.g. exploiting macrophage phagocytosis). Minority variants are highly concentrated around the gene portion encoding for the Spike cleavage site, with a higher incidence in the URT; four mutations are highly recurring among samples and were found associated with the URT. Interestingly, 55.8% of minority variants detected in this locus were T>G and G>T transversions. Results from this study evidenced the presence of selective pressure and suggest that an evolutionary process is still ongoing in one of the crucial sites of spike protein associated with the spillover to humans.

SARS-CoV-2 Molecular Epidemiology Can Be Enhanced by Occupational Health: The Experience of Monitoring Variants of Concern in Workplaces in Rio de Janeiro, Brazil.

The emergence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to extra caution in workplaces to avoid the coronavirus disease 2019 (COVID-19). In the occupational environment, SARS-CoV-2 testing is a powerful approach in providing valuable information to detect, monitor, and mitigate the spread of the virus and preserve productivity. Here a centralized Occupational Health Center provided molecular diagnosis and genomic sequences for companies and industries in Rio de Janeiro, Brazil. From May to August 2021, around 20% of the SARS-CoV-2 positive nasopharyngeal swabs from routinely tested workers were sequenced and reproduced the replacement of Gamma with Delta variant observed in regular surveillance programs. Moreover, as a proof-of-concept on the sensibility of the occupational health genomic surveillance program described here, it was also found: i) the primo-identification of B.1.139 and A.2.5 viral genomes in Brazil and ii) an improved dating of Delta VoC evolution, by identifying earlier cases associated with AY-related genomes. We interpret that SARS-CoV-2 molecular testing of workers, independent of symptom presentation, provides an earlier opportunity to identify variants. Thus, considering the continuous monitoring of SARS-CoV-2 in workplaces, positive samples from occupation health programs should be regarded as essential to improve the knowledge on virus genetic diversity and VoC emergence.

Rapid Replacement of SARS-CoV-2 Variants by Delta and Subsequent Arrival of Omicron, Uganda, 2021.

Genomic surveillance in Uganda showed rapid replacement of severe acute respiratory syndrome coronavirus 2 over time by variants, dominated by Delta. However, detection of the more transmissible Omicron variant among travelers and increasing community transmission highlight the need for near-real-time genomic surveillance and adherence to infection control measures to prevent future pandemic waves.

Rapid SARS-CoV-2 Adaptation to Available Cellular Proteases.

Recent emergence of SARS-CoV-1 variants demonstrates the potential of this virus for targeted evolution, despite its overall genomic stability. Here we show the dynamics and the mechanisms behind the rapid adaptation of SARS-CoV-2 to growth in Vero E6 cells. The selective advantage for growth in Vero E6 cells is due to increased cleavage efficiency by cathepsins at the mutated S1/S2 site. S1/S2 site also constitutes a heparan sulfate (HS) binding motif that influenced virus growth in Vero E6 cells, but HS antagonist did not inhibit virus adaptation in these cells. The entry of Vero E6-adapted virus into human cells is defective because the mutated spike variants are poorly processed by furin or TMPRSS2. Minor subpopulation that lack the furin cleavage motif in the spike protein rapidly become dominant upon passaging through Vero E6 cells, but wild type sequences are maintained at low percentage in the virus swarm and mediate a rapid reverse adaptation if the virus is passaged again on TMPRSS2+ human cells. Our data show that the spike protein of SARS-CoV-2 can rapidly adapt itself to available proteases and argue for deep sequence surveillance to identify the emergence of novel variants. IMPORTANCE Recently emerging SARS-CoV-2 variants B.1.1.7 (alpha variant), B.1.617.2 (delta variant), and B.1.1.529 (omicron variant) harbor spike mutations and have been linked to increased virus pathogenesis. The emergence of these novel variants highlights coronavirus adaptation and evolution potential, despite the stable consensus genotype of clinical isolates. We show that subdominant variants maintained in the virus population enable the virus to rapidly adapt to selection pressure. Although these adaptations lead to genotype change, the change is not absolute and genomes with original genotype are maintained in the virus swarm. Thus, our results imply that the relative stability of SARS-CoV-2 in numerous independent clinical isolates belies its potential for rapid adaptation to new conditions.

Multiorgan and Vascular Tropism of SARS-CoV-2.

Although the respiratory tract is the main target of SARS-CoV-2, other tissues and organs are permissive to the infection. In this report, we investigated this wide-spectrum tropism by studying the SARS-CoV-2 genetic intra-host variability in multiple tissues. The virological and histological investigation of multiple specimens from a post-mortem COVID-19 patient was performed. SARS-CoV-2 genome was detected in several tissues, including the lower respiratory system, cardio-vascular biopsies, stomach, pancreas, adrenal gland, mediastinal ganglion and testicles. Subgenomic RNA transcripts were also detected, in favor of an active viral replication, especially in testicles. Ultra-deep sequencing allowed us to highlight several SARS-CoV-2 mutations according to tissue distribution. More specifically, mutations of the spike protein, i.e., V341A (18.3%), E654 (44%) and H655R (30.8%), were detected in the inferior vena cava. SARS-CoV-2 variability can contribute to heterogeneous distributions of viral quasispecies, which may affect the COVID-19 pathogeny.

SARS-CoV-2 variant detection with ADSSpike.

The SARS-CoV-2 coronavirus pandemic has been an unprecedented challenge to global pandemic response and preparedness. With the continuous appearance of new SARS-CoV-2 variants, it is imperative to implement tools for genomic surveillance and diagnosis in order to decrease viral transmission and prevalence. The ADSSpike workflow was developed with the goal of identifying signature SNPs from the S gene associated with SARS-CoV-2 variants through amplicon deep sequencing. Seventy-two samples were sequenced, and 30 mutations were identified. Among those, signature SNPs were linked to 2 Zeta-VOI (P.2) samples and one to the Alpha-VOC (B.1.17). An average depth of 700 reads was found to properlycorrectly identify all SNPs and deletions pertinent to SARS-CoV-2 mutants. ADSSpike is the first workflow to provide a practical, cost-effective, and scalable solution to diagnose SARS-CoV-2 VOC/VOI in the clinical laboratory, adding a valuable tool to public health measures to fight the COVID-19 pandemic for approximately $41.85 USD/reaction.